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ABSTRACT 

, The BL Lac object 3C 66A was the target of an extensive multiwavelength monitoring cam- 

paign from July 2003 through April 2004 (with a core campaign from Sept. - Dec. 2003), 
involving observations throughout the electromagnetic spectrum. Radio, infrared, and optical 
observations were carried out by the WEBT-ENIGMA collaboration. At higher energies, 3C 66A 
was observed in X-rays (RXTE), and at very-high-energy (VHE) 7 -rays (STACEE, VERITAS). 
In addition, the source has been observed with the VLBA at 9 epochs throughout the period 
J> , September 2003 - December 2004, including 3 epochs contemporaneous with the core campaign. 

■ A gradual brightening of the source over the course of the campaign was observed at all optical 

frequencies, culminating in a very bright maximum around Feb. 18, 2004. The WEBT campaign 
revealed microvariability with flux changes of ^ 5 % on time scales as short as ~ 2 hr. The 
source was in a relatively bright state, with several bright flares on time scales of several days. 
The spectral energy distribution (SED) indicates a vF v peak in the optical regime. A weak trend 
of optical spectral hysteresis with a trend of spectral softening throughout both the rising and 
decaying phases, has been found. On longer time scales, there appears to be a weak indication of 
a positive hardness-intensity correlation for low optical fluxes, which does not persist at higher 
flux levels. 

The 3-10 keV X-ray flux of 3C 66A during the core campaign was historically high and its 
spectrum very soft, indicating that the low- frequency component of the broadband SED extends 
beyond ~ 10 keV. No significant X-ray flux and/or spectral variability was detected. STACEE 
and Whipple observations provided upper flux limits at > 150 GeV and > 390 GeV, respectively. 

The 22 and 43 GHz data from the 3 VLBA epochs made between September 2003 and 
January 2004 indicate a rather smooth jet with only very moderate internal structure. Evidence 
for superluminal motion (8.5 ± 5.6 ft. -1 c) was found in only one out of 6 components, while the 
apparent velocities of all other components are consistent with 0. The radial radio brightness 
profile suggests a magnetic field decay cx r -1 and, thus, a predominantly perpendicular magnetic 
field orientation. 
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1. Introduction 

BL Lac objects and flat-spectrum radio quasars 
(FSRQs) are active galactic nuclei (AGNs) com- 
monly unified in the class of blazars. They exhibit 
some of the most violent high-energy phenomena 
observed in AGNs to date. They are characterized 
by non-thermal continuum spectra, a high degree 
of linear polarization in the optical, rapid variabil- 
ity at all wavelengths, radio jets with individual 
components often exhibiting apparent superlumi- 
nal motion, and — at least episodically — a signif- 
icant portion of the bolometric flux emitted in 7- 
rays. 46 blazars have been detected and identified 
with high confidence in high energy (> 100 MeV) 
gamma-rays by the EGRET instrument on board 
the Compton Gamma-Ray Observatory (Hartman 
et al. 1999; Mattox, Hartman, & Reimer 2001). 
The properties of BL Lac objects and blazar-type 
FSRQs are generally very similar, except that 
BL Lac objects usually show only weak emission 
or absorption lines (with equivalent width in the 
rest-frame of the host galaxy of < 5 A) , if any. In 
3C 66A (= 0219+428 = NRAO 102 = 4C 42.07), 
a weak Mg II emission line has been detected by 
Miller, French, & Hawley (1978). This led to the 
determination of its rcdshift at z = 0.444, which 
was later confirmed by the detection of a weak Lya 
line in the IUE spectrum of 3C 66A (Lanzetta, 
Turnshek, & Sandoval 1993). However, as recently 
pointed out by Bramel et al. (2005), these redshift 
determinations are actually still quite uncertain. 
In this paper, we do base our analysis on a redshift 
value of z = 0.444, but remind the reader that 
some results of the physical interpretation should 
be considered as tentative pending a more solid 
rcdshift determination. 

In the framework of relativistic jet models, the 
low- frequency (radio - optical/UV) emission from 
blazars is interpreted as synchrotron emission from 
nonthermal electrons in a relativistic jet. The 
high-frequency (X-ray - 7-ray) emission could ei- 
ther be produced via Compton upscattering of 
low frequency radiation by the same electrons re- 
sponsible for the synchrotron emission (lcptonic 
jet models; for a recent review see, e.g., Bottcher 
2002), or due to hadronic processes initiated by 
relativistic protons co-accelerated with the elec- 
trons (hadronic models, for a recent discussion see, 
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e.g., Miicke & Protheroe 2001; Miicke et al. 2003). 

To date, 6 blazars have been detected at very 
high energies (> 300 GeV) with ground-based air 
Cerenkov detectors (Punch et al. 1992; Quinn et 
al. 1996; Catanese et al. 1998; Chadwick et al. 
1999; Aharonian et al. 2002; Horan et al. 2002; 
Holder et al. 2003). All of these belong to the 
sub-class of high-frequency peaked BL Lac ob- 
jects (HBLs). The field of extragalactic GeV - 
TeV astronomy is currently one of the most excit- 
ing research areas in astrophysics, as the steadily 
improving flux sensitivities of the new generation 
of air Cerenkov telescope arrays and their de- 
creasing energy thresholds (for a recent review 
see, e.g., Weekes et al. 2002), provides a grow- 
ing potential to extend their extragalactic-source 
list towards intermediate and even low-frequency 
peaked BL Lac objects (LBLs) with lower vF v 
peak frequencies in their broadband spectral en- 
ergy distributions (SEDs). Detection of such ob- 
jects at energies ~ 40 - 100 GeV might provide 
an opportunity to probe the intrinsic high-energy 
cutoff of their SEDs since at those energies, 77 
absorption due to the intergalactic infrared back- 
ground is still expected to be small (e.g., de Jager 
& Stecker 2002). 

3C 66A has been suggested as a promising 
candidate for detection by the new generation 
of atmospheric Cerenkov telescope facilities like 
STACEE or VERITAS (e.g. Costamante & Ghis- 
ellini 2002). Neshpor et al. (1998, 2000) have ac- 
tually reported multiple detections of the source 
with the GT-48 telescope of the Crimean Astro- 
physical Observatory, but those could not be con- 
firmed by any other group so far (see, e.g. Horan 
et al. 2004). 

3C 66A is classified as a low-frequency peaked 
BL Lac object (LBL), a class also commonly re- 
ferred to as radio selected BL Lac objects. Its 
low-frequency spectral component typically peaks 
at IR - UV wavelengths, while the high-frequency 
component seems to peak in the multi-MeV - GeV 
energy range. Since its optical identification by 
Wills & Wills (1974), 3C 66A has been the target 
of many radio, IR, optical, X-ray, and 7-ray ob- 
servations in the past, although it is not as regu- 
larly monitored at radio frequencies as many other 
blazars due to problems with source confusion 
with the nearby radio galaxy 3C 66B (6 '.5 from 
3C 66 A), in particular at lower (4.8 and 8 GHz) 
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frequencies (Aller, Aller, & Hughes 1994; Takalo 
et al. 1996). 

The long-term variability of 3C 66A at near- 
infrared (J, H, and K bands) and optical (U, 
B, V, R, I) wavelengths has recently been com- 
piled and analyzed by Fan & Lin (1999) and Fan 
& Lin (2000), respectively. Variability at those 
wavelengths is typically characterized by varia- 
tions over < 1.5 mag on time scales ranging from 
~ 1 week to several years. A positive correla- 
tion between the B - R color (spectral hardness) 
and the R magnitude has been found by Vagnetti 
et al. (2003). An intensive monitoring effort by 
Takalo et al. (1996) revealed evidence for rapid 
microvariability, including a decline ~ 0.2 mag 
within ~ 6 hr. Microvariability had previously 
been detected in 3C 66A by Carini & Miller (1991) 
and De Diego et al. (1997), while other, similar 
efforts did not reveal such evidence (e.g. Miller 
& McGimsey 1978; Takalo et al. 1992; Xie et al. 
1992). Lainela et al. (1999) also report on a 65- 
day periodicity of the source in its optically bright 
state, which has so far not been confirmed in any 
other analysis. 

3C 66A is generally observed as a point source, 
with no indication of the host galaxy. The 
host galaxy of 3C 66A was marginally resolved 
by Wurtz, Stocke, & Yee (1996). They found 
i?Gunn = 19.0 mag for the host galaxy; the Hubble 
type could not be determined. 

In X-rays, the source has been previously de- 
tected by EXOSAT (Sambruna et al. 1994), EIN- 
STEIN (Worrall & Wilkes 1990), ROSAT (Fos- 
sati et al. 1998), BeppoSAX (Perri et al. 2003), 
and XMM-Newton (Croston et al. 2003). It shows 
large-amplitude soft-X-ray variability among these 
various epochs of observation, with flux levels at 
1 keV ranging from ~ 0.4 /iJy to ~ 5 /iJy and 
generally steep (energy index a > 1) soft X-ray 
spectra below 1 keV. 3C 66A has also been de- 
tected in > 100 MeV 7-rays by EGRET on sev- 
eral occasions, with flux levels up to F>iooMeV = 
(25.3 ± 5.8) x 10~ 8 photons cm~ 2 s" 1 (Hartman 
et al. 1999). 

Superluminal motion of individual radio com- 
ponents of the jet has been detected by Jorstad et 
al. (2001). While the identification of radio knots 
across different observing epochs is not unique, 
Jorstad et al. (2001) favor an interpretation imply- 
ing superluminal motions of up to /3 ap p ~ 19 h^ 1 « 



27. This would imply a lower limit on the bulk 
Lorentz factor of the radio emitting regions of 
T > 27. However, theoretical modeling of the 
non-simultaneous SED of 3C 66A (Ghiscllini et al. 
1998) suggests a bulk Lorentz factor of the emit- 
ting region close to the core — where the 7-ray 
emission is commonly believed to be produced — 
of T ~ 14, more typical of the values obtained for 
other blazars as well. 

In spite of the considerable amount of observa- 
tional effort spent on 3C 66A (see, e.g. Takalo et al. 
1996, for an intensive optical monitoring campaign 
on this source), its multiwavelength SED and cor- 
related broadband spectral variability behaviour 
are still surprisingly poorly understood, given its 
possible VHE 7-ray source candidacy. The object 
has never been studied in a dedicated multiwave- 
length campaign during the life time of EGRET. 
There have been few attempts of coordinated mul- 
tiwavelength observations. For example, Worrall 
et al. (1984) present quasi-simultaneous radio, IR, 
optical, and UV observations of 3C 66A in 1983, 
but observations at different wavelength bands 
were still separated by up to ~ 2 weeks, and no 
simultaneous higher-energy data were available. 
This is clearly inadequate to seriously constrain re- 
alistic, physical emission models, given the estab- 
lished large-amplitude variability on similar time 
scales and the importance of the high-energy emis- 
sion in the broadband SED of the source. 

For this reason, we have organized an inten- 
sive multiwavelength campaign to observe 3C 66A 
from July 2003 through April 2004, focusing on 
a core campaign from Sept. - Dec. 2003. In 
§2, we describe the observations and data analysis 
and present light curves in the various frequency 
bands. Spectral variability patterns are discussed 
in §3, and the results of our search for inter-band 
cross-correlations and time lags are presented in 
§4. Simultaneous broadband spectral energy dis- 
tributions (SEDs) of 3C 66A at various optical 
brightness levels are presented in §5. In §6 we 
use our results to derive estimates of physical pa- 
rameters, independent of the details of any spe- 
cific model. We summarize in §7. In a companion 
paper (Joshi & Bottcher, 2005, in preparation), 
we will use time-dependent leptonic models to fit 
the spectra and variability patterns found in this 
campaign, and make specific predictions concern- 
ing potentially observable X-ray spectral variabil- 
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ity patterns and 7-ray emission. 

Throughout this paper, we refer to a as the 
energy spectral index, F v [Jy] oc v~ a . A cos- 
mology with fl m = 0.3, Oa = 0.7, and Hq = 
70 km s _1 Mpc -1 is used. In this cosmology, and 
using the redshift of z = 0.444, the luminosity dis- 
tance of 3C 66A is d L = 2.46 Gpc. 

2. Observations, data reduction, and light 
curves 

3C 66A was observed in a coordinated mul- 
tiwavelength campaign at radio, near-IR, opti- 
cal (by the WEBT- ENIGMA collaboration), X- 
ray, and VHE 7-ray energies during a core cam- 
paign period of September 2003 - December 2003. 
The object is being continuously monitored at 
radio and optical wavelength by several ongoing 
projects, and we present data in those wavelengths 
regimes starting in July 2003. During the core 
campaign, wc found that the source was grad- 
ually brightening, indicating increasing activity 
throughout the campaign period. For this reason, 
the WEBT campaign was extended, and we con- 
tinued to achieve good time coverage until early 
March 2004, and individual observatories still con- 
tributed data through April 2004. The overall 
timeline of the campaign, along with the mea- 
sured long-term light curves at radio, infrared, op- 
tical, and X-ray frequencies is illustrated in Fig. 1. 
Table 1 lists all participating observatories which 
contributed data to this campaign. In this section, 
we will describe the individual observations in the 
various frequency ranges and outline the data re- 
duction and analysis. 

2.1. Radio observations 

At radio frequencies, the object was monitored 
using the University of Michigan Radio Astron- 
omy Observatory (UMRAO) 26 m telescope, at 
4.8, 8, and 14.5 GHz, the 14 m Metsahovi Radio 
Telescope of the Helsinki University of Technol- 
ogy, at 22 and 37 GHz, and the 576-m ring tele- 
scope (RATAN-600) of the Russian Academy of 
Sciences, at 2.3, 3.9, 7.7, 11, and 22 GHz. 

At the UMRAO, the source was monitored in 
the course of the on-going long-term blazar moni- 
toring program. The data were analyzed following 
the standard procedure described in Aller et al. 
(1985). As mentioned above, the sampling at the 



lower frequencies (4.5 and 8 GHz) was rather poor 
(about once every 1-2 weeks) and some individ- 
ual errors were rather large due to source confusion 
problems with 3C 66B. At 14.5 GHz, a slightly 
better sampling of ~ 1 - 2 observations per week, 
at least until the end of 2003, was achieved with 
relative flux errors of typically a few percent. The 
resulting 14.5 GHz light curve is included in Fig. 1. 
It seems to indicate generally moderate variabil- 
ity (AF/F < 25 %) on time scales of > 1 week, 
though a discrete auto-correlation analysis (Edel- 
son & Krolik 1988) does not reveal any significant 
structure, primarily due to the insufficient sam- 
pling. 

At 22 and 37 GHz, 3C 66A has been monitored 
using the 14 m radio telescope of the Metsahovi 
Radio Observatory of the Helsinki University of 
Technology. The data have been reduced with 
the standard procedure described in Terasranta et 
al. (1998). The 22 and 37 GHz light curves re- 
veal moderate-amplitude (AF/F < 30 %), erratic 
variability which is clearly under-sampled by the 
available data. The discrete auto-correlation func- 
tions indicate a time scale of a few days for this 
short-term variability. 

We point out that the observed erratic variabil- 
ity may, at least in part, be a consequence of in- 
terstellar scintillation. At the Galactic coordinates 
of the source, the transition frequency, where the 
interstellar scattering strength (a measure of the 
phase change induced by interstellar scattering) 
becomes unity, is - 7 GHz (Walker 1998). At 
higher frequencies, as considered here, scattering 
will occur in the weak scattering regime, with only 
small phase changes through interstellar scatter- 
ing. In this regime, we find the fractional rms vari- 
ability amplitude for a point source due to inter- 
stellar scintillation as 0.36 and 0.09 at a frequency 
of 14.5 and 37 GHz, respectively, and the respec- 
tive variability time scales due to interstellar scin- 
tillation are 1.4 and 0.9 hours. Consequently, in 
particular at lower frequencies (< 20 GHz), a sub- 
stantial fraction of the observed variability may 
well be due to interstellar scintillation. 

The RATAN-600 has monitored 3C 66A, per- 
forming 17 observations from mid-1997 through 
Oct. 2003. The last two of these observations 
(Oct. 11 and 14) coincided with our core cam- 
paign, and provided additional frequency cover- 
age at 2.3, 3.9, 7.7, 11, and 22 GHz. Those data 
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Table 1 

List of observatories that contributed data to this campaign 



Observatory 


Specifications 


frequency / niters / energy range 


N aba 


Radio Observatories 


UMRAO, Michigan, USA 


26 m 




4.8, 8, 14.5 GHz 


46 


Mctsahovi, Finland 


14 m 




22, 37 GHz 


103 


RATAN-600, Russia 


576m (ring) 


2.3 


, 3.9, 7.7, 11, 22 GHz 


17 


VLBA 


10x25 m 


2.3, 


5, 8.4, 22, 43, 86 GHz 


9 


Infrared Observatories 


Campo Impcratore, Italy 


1.1 m 




J, H, K 


171 


Mt. Abu, India (MIRO) 


1.2 m (NICMOS-3) 




J, H, K' 


79 


Roque (NOT), Canary Islands 


2.56 m 




J, H, K 


15 


Optical Observatories 


Abastumani, Georgia (FSU) 


70 cm 




R 


210 


Armenzano, Italy 


40 cm 




B, V, R, I 


315 


Bell Obs., Kentucky, USA 


60 cm 




R 


12 


Boltwood, Canada 


40 cm 




B, V, R, I 


402 


Catania, Italy 


91 cm 




U, B, V 


835 


Crimean Astr. Obs., Ukraine 


70 cm ST-7 




B, V, R, I 


85 


Heidelberg, Germany 


70 cm 




B, R, I 


8 


Kitt Peak (MDM), Arizona, USA 


130 cm 




B, V, R, I 


147 


Michael Adrian Obs., Germany 


120 cm 




R 


30 


Mt. Lemmon, Korea 


100 cm 




B, V, R, I 


399 


Mt. Maidanak, Uzbekistan 


150 cm AZT-22 




B, R 


1208 


Nyrola, Finland 


40 cm SCT 




R 


159 


Osaka Kyoiku, Japan 


51 cm 




R 


1167 


Perugia, Italy 


40 cm 




V, R, I 


140 


Roque (KVA), Canary Islands 


35 cm 




B, V, R 


653 


Roque (NOT), Canary Islands 


256 cm 




U, B, V, R, I 


52 


Sabadcll, Spain 


50 cm 




B, R 


4 


San Pedro Martfr, Mexico 


150 cm 




B, V, R, I 


185 


Shanghai, China 


156 cm 




V, R 


36 


Skinakas, Crete 


130 cm 




B, V, R, I 


156 


Sobacksan, Korea 


61 cm 




B, V, R, I 


133 


St. Louis, Missouri, USA 


35 cm 




B, R 


16 


Torino, Italy 


105 cm REOSC 




B, V, R, I 


227 


Tuorla, Finland 


103 cm 




B, V, R 


1032 


X-R,ay Observatory 


RXTE 


PCA 




3-25 keV 


26 


7-R.ay Observatories 


STACEE 


Solar-Tower Cerenkov Array 




> 100 GeV 


85 


VERITAS 


Whipple 10 m 




> 390 GeV 


31 
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were analyzed as described in detail in Kovalev et 
al. (1999), and the time average of the resulting 
fluxes from the two observations in Oct. 2003 are 
included in the spectral energy distribution shown 
in fig. 15. 

A total of 9 observing epochs using the VLBA 2 
have been approved to accompany our multiwave- 
length campaign on 3C 66A. In this paper, we de- 
scribe details of the data analysis and results of 
the first (2003.78; VLBA observing code BS133A), 
the second (2003.83; BS133B), and the fourth 
(2004.08; BS133E) epoch, concentrating on the re- 
sults at 22 and 43 GHz. The third epoch (2003.96) 
suffered from radio frequency interference, and its 
reduction is in progress. Results of the complete 
set of all 9 VLBA observing epochs, including all 
6 frequencies (2.3, 5, 8.4, 22, 43, and 86 GHz) to- 
gether with polarization data, will be presented in 
a separate paper (Savolainen et al., 2005, in prepa- 
ration). 

The data were correlated on the VLBA correla- 
tor and were postprocessed with the NRAO Astro- 
nomical Image Processing System, AIPS (Greisen 
1988), and the Caltech DIFMAP package (Shep- 
herd 1997). Standard methods for VLBI data 
reduction and imaging were used. A priori am- 
plitude calibration was performed using measured 
system temperatures and gain curves in the AIPS 
task APCAL. At this point, a correction for atmo- 
spheric opacity was also applied. After removal of 
the parallactic angle phase, single-band delay and 
phase offsets were calculated manually by fringe- 
fitting a short scan of data of a bright calibrator 
source (0420-014). We did manual phase calibra- 
tion instead of using pulsc-cal tones, because there 
were unexpected jumps in the phases of the pulse- 
cal during the observations. Global fringe-fitting 
was performed with the AIPS task FRING using 
Los Alamos (LA) as a reference antenna. The 
AIPS procedure CRSFRING was used to remove 
the delay difference between right- and left-hand 
systems (this is needed for polarization calibra- 
tion). A bandpass correction was determined and 
applied before averaging across the channels, after 
which the data were imported into DIFMAP. 



2 The VLBA is a facility of the National Radio Astron- 
omy Observatory (NRAO). The NRAO is a facility of the 
National Science Foundation, operated under cooperative 
agreement by Associated Universities, Inc. 
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Fig. 1. — Timeline of the broadband campaign 
on 3C 66 A in 2003/2004. The grey vertical lines 
indicate the epochs of the VLBA observations. 



3C 88A at 33 GHz — gaussian models 




Relative R.A. (marcsec) 



Fig. 2.— 22-GHz VLBA maps of 3C 66A dur- 
ing the three epochs in 2003/2004 contemporane- 
ous with our campaign. Overlaid are the fit re- 
sults from decomposing the structure into a sum 
of Gaussian components. 



7 



In DIFMAP, the data were first phase self- 
calibrated using a point source model, and then 
averaged in time. We performed data editing in a 
station-based manner, and ran several iterations 
of CLEAN and phase self-calibration in Stokes 
I. After a reasonable fit to the closure phases 
was obtained, we also performed amplitude self- 
calibration, first with a solution interval corre- 
sponding to the whole observation length. The 
solution interval was then gradually shortened as 
the model improved by further CLEANing. Final 
images were produced using the Perl library FIT- 
Splot 3 . For final images, normal weighting of the 
uv-data was used in order to reveal the extended, 
low surface brightness emission. 

We have checked the amplitude calibration by 
comparing the extrapolated zero baseline flux den- 
sity of the compact source 0420-014 to the single- 
dish flux measurements at Metsahovi. The fluxes 
are comparable: at 22 GHz the VLBI flux is on 
average 5 % lower than the single-dish flux, and 
the 43 GHz flux is on average 8 % lower than 
the nearly simultaneous single-dish measurements 
at 37 GHz. The small amount of missing flux is 
probably due to some extended emission resolved 
out in the high frequency VLBI images. The inte- 
grated VLBI flux of 3C 66A is about 30 % smaller 
than the corresponding single-dish value both at 
22 GHz and at 43 GHz, but since the source has 
notable kpc-scale structure evident in the VLA im- 
ages (Price et al. 1993; Taylor et al. 1996), it is 
most likely that the missing flux comes from this 
kpc-scale jet which is resolved out in our VLBI 
images. Thus, we conclude that the accuracy of 
the amplitude calibration is better than 10 % at 
both frequencies. 

In order to estimate the parameters of the emis- 
sion regions in the jet, we did model fitting to the 
self-calibrated visibilities in DIFMAP. The data 
were fitted with circular Gaussian model compo- 
nents and we sought to obtain the best possible fit 
to the visibilities and to the closure phases. Sev- 
eral starting points were tried in order to avoid a 
local minimum fit. The results of this fitting pro- 
cedure are included in Figs. 2 and 3. In all the 
model-fits, the core of the VLBI jet is the north- 
ernmost and also the brightest component. Be- 
yond the core, we have divided the jet into three 

http : //personal . denison . edu/~homand/ 



regions named A, B and C. Closest to the core is 
the region C, which consists of three components 
with decreasing surface brightness as a function 
of the distance from the core. The region B, also 
made of three components, shows a clear bending 
of the jet together with re-brightening. The ob- 
served increase in the flux at this point could plau- 
sibly be due to increased Doppler boosting caused 
by the bending of the jet towards our line of sight. 
The region A shows weak extended emission at 8 
mas from the core at 22 GHz. This emission is 
more pronounced at lower frequencies. 

Fig. 4 shows the separation of the model com- 
ponents from the core, at each epoch. The com- 
ponents A - C are weak compared to the core 
(see Table 3), which renders the estimation of er- 
rors of the model parameters using programs like 
Difwrap (Lovell 2000) problematic. Thus, in Fig. 
4, we have assumed the following uncertainties in 
the component position accuracy: 1) for bright 
features with SNR > 50, the positions are ac- 
curate to 20 % of the projection of the ellipti- 
cal beam FWHM onto a line joining the center 
of the core to the center of the components; 2) 
for weaker components with SNR < 50, uncer- 
tainties are 50 % of the beam size or 50 % of the 
size of the knot — whichever is larger. These es- 
timates are based on the experience with other 
sources (Savolainen et al. 2005, in preparation; 
Wiik et al. 2005, in preparation), and they follow 
the results by Jorstad et al. (2005), where uncer- 
tainties of model-fit parameters for a large number 
of observations are estimated. Positional errors 
of similar size are also reported by Homan et al. 
(2001), who have estimated the uncertainties from 
the variance of the component position about the 
best-fit polynomial. 

There is a ~ 0.1 mas shift of the components CI 
- C3 between the two frequencies, with the 22 GHz 
model components appearing further downstream. 
This cannot be explained by an opacity effect, 
since in that case the components would have 
shifted in the opposite direction: at higher fre- 
quency we would expect to see emission from the 
region closer to the apex of the jet than at lower 
frequency, if the opacity effect is significant. This 
obviously is not the case in Fig. 4, and more- 
over, for the components B2 and B3 the positions 
at 22 and 43 GHz arc coincident. However, there 
is a simple explanation for the observed shift, if 
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Table 2 

Core brightness temperatures and corresponding lower limits, calculated from the 
measured 43 ghz vlba fluxes and component sizes 



Epoch T§°™ [K] [K] 

2003.78 2.7 x 10 11 1.1 x 10 11 

2003.83 2.0 x 10 11 1.0 x 10 11 

2004.08 6.6 x 10 11 1.4 x 10 11 



Table 3 

vlba component fluxes, core distances (rad.), component position angles (pa), and 

circular component diameters (size) 



Epoch Freq. [GHz] Comp. ID Flux [mjy] Rad. [mas] PA [deg] Size [mas] 



2003.78 


22 GHZ 


Core 


387 


0.0 


0.0 


0.047 






C3 


61 


0.295 


-172.9 


0.136 






C2 


36 


0.658 


-168.2 


0.164 






CI 


19 


1.096 


-167.9 


0.309 






B3 


21 


1.803 


-166.1 


0.278 






B2 


40 


2.287 


-174.6 


0.331 






Bl 


19 


3.086 


-176.2 


0.970 






A 


37 


8.439 


172.8 


3.333 




43 GHz 


Core 


355 


0.0 


0.0 


0.029 






C3 


63 


0.183 


-174.0 


0.089 






C2 


43 


0.534 


-171.1 


0.218 






CI 


24 


0.940 


-165.4 


0.304 






B3 


18 


1.849 


-167.3 


0.284 






B2 


30 


2.279 


-173.8 


0.329 






Bl 


9 


2.821 


-173.8 


0.349 


2003.83 


22 GHz 


Core 


375 


0.0 


0.0 


0.032 






C3 


59 


0.282 


-172.4 


0.118 






C2 


38 


0.685 


-170.2 


0.193 






CI 


18 


1.167 


-167.3 


0.298 






B3 


23 


1.892 


-167.7 


0.274 






B2 


39 


2.320 


-174.8 


0.358 






Bl 


13 


3.051 


-177.6 


0.836 






A 


41 


8.546 


169.7 


3.836 




43 GHz 


Core 


369 


0.0 


0.0 


0.034 






C3 


76 


0.187 


-176.3 


0.116 






C2 


34 


0.584 


-172.1 


0.210 






CI 


26 


1.050 


-166.0 


0.383 






B3 


26 


1.947 


-168.6 


0.385 






B2 


20 


2.328 


-175.0 


0.213 






Bl 


16 


2.829 


-174.1 


0.795 


2004.08 


22 GHz 


Core 


363 


0.0 


0.0 


0.046 






C3 


64 


0.289 


-172.4 


0.132 






C2 


32 


0.712 


-172.0 


0.235 






CI 


12 


1.241 


-167.2 


0.311 






B3 


22 


1.860 


-165.7 


0.296 






B2 


45 


2.333 


-174.2 


0.386 






Bl 


18 


3.285 


-178.5 


1.754 






A 


38 


8.348 


171.9 


3.702 




43 GHz 


Core 


372 


0.0 


0.0 


0.019 






C3 


59 


0.207 


-171.2 


0.078 






C2 


30 


0.568 


-171.6 


0.207 






CI 


13 


1.109 


-166.2 


0.252 






B3 


15 


1.816 


-166.3 


0.320 






B2 


40 


2.289 


-173.7 


0.411 






XI 


4 


3.375 


-171.8 


0.417 
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the brightness profile along the jet in region C is 
smooth. Namely, the flux in the outer part of 
the region C decreases more steeply at 43 GHz 
than at 22 GHz (see Figs. 2 and 3). As the re- 
gion C is modelled by the same number of compo- 
nents at both frequencies, the model-fit procedure 
shifts the 22 GHz components downstream rela- 
tive to 43 GHz components in order to represent 
the power-law. 

Fig. 4 shows that, except for one component, 
all the components revealed by the analysis of the 
three epochs considered here are consistent with 
zero proper motion. Although a monitoring effort 
over a longer time scale will be necessary to settle 
this issue, our analysis shows slower component 
motion than what is presented by Jorstad et al. 
(2001), where fast superluminal motions of several 
components, with speeds ranging from 9 h^ 1 c to 
19 ft -1 c were found. Note, however, that the data 
used by Jorstad et al. (2001) covered a much longer 
time frame and refers to different components and 
a rather different state of activity of 3C 66A. More- 
over, there are also similarities between the data 
sets: Jorstad et al. (2001) reported a stationary 
component (their component "C") at distance of 
^0.5 mas from the core, which could correspond 
to our component C2. New 7 mm VLB A monitor- 
ing data presented in a recent paper by Jorstad et 
al. (2005) suggest that there are two kinds of com- 
ponents in 3C 66A: fast and very weak components 
with apparent speeds > 20 c, and stronger com- 
ponents with moderate velocities of 1.5 - 5 c. The 
components CI - C3 in our analysis could be qual- 
itatively similar to the latter. More solid conclu- 
sions concerning the question of superluminal mo- 
tions of individual components might be possible 
after the final analysis of all 9 epochs of the entire 
VLBA monitoring program of 3C66A proposed in 
connection with this campaign (Savolainen et al., 
2005, in preparation). Based on the analysis of the 
three epochs contemporaneous with the other ob- 
servations of this campaign, one component, CI, 
shows superluminal motion of (8.5 ± 5.6) h^ 1 c = 
(12.1 ±8.0) c for h = 0.7. 

The essentially zero proper motion of most of 
the components together with the observed trend 
of the brightness temperature vs. distance along 
the jet (see §6.2) suggest that in the case of 3C 66A 
the VLBA model-fit components might not cor- 
respond to actual 'knots' in the jet, but a rather 



3C 66A at 43 GHz - gaussian i 




Fig. 3.— 43-GHz VLBA maps of 3C 66A dur- 
ing the three epochs in 2003/2004 contemporane- 
ous with our campaign. Overlaid are the fit re- 
sults from decomposing the structure into a sum 
of Gaussian components. 
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Fig. 4. — Separation from the core of the Gaus- 
sian fit components shown in Figs. 2 and 3, as a 
function of time. Filled symbols refer to 22 GHz, 
open symbols to 43 GHz. The solid (dashed) lines 
show linear least-squares fits to the positions of 
22 GHz (43 GHz) components. The velocity re- 
sults are the average of the results from the two 
frequencies. All components except CI are consis- 
tent with being stationary. 
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smooth flow with simple brightness profile (such as 
a power law) might better describe the jet during 
this campaign — at least in region C. More de- 
tailed parameter estimates will be extracted from 
these results in §6.2. 

2.2. Infrared observations 

In the context of the extensive WEBT cam- 
paign (see §2.3), 3C 66A was also observed at 
near-infrared wavelengths in the J, H, and K/K' 
bands at three observatories: The Campo Impera- 
tore 1.1-m telescope of the Infrared Observatory of 
the Astronomical Observatory of Rome, Italy, the 
1.2-m telescope (using the NICMOS-3 HgCdTe IR 
array with 256 x 256 pixels, with each pixel cor- 
responding to 0.96" on the sky) at the Mt. Abu 
Infrared Observatory (MIRO) at Mt. Abu, In- 
dia, and at the 2.56-m Nordic Optical Telescope 
(NOT) on Roque de los Muchachos on the Canary 
Island of La Palma. The primary data were an- 
alyzed using the same standard technique as the 
optical data (see §2.3), including flat-field subtrac- 
tion, extraction of instrumental magnitudes, cali- 
bration against comparison stars to obtain stan- 
dard magnitudes, and de-reddening. The sam- 
pling was generally not dense enough to allow an 
improvement of the data quality by re-binning. 
Unfortunately, the three observatories did not per- 
form any observations on the same day, so that the 
cross-calibration between different instruments is 
problematic. We have therefore opted not to cor- 
rect for possible systematic offsets beyond the cal- 
ibration to standard magnitudes of well-calibrated 
comparison stars (see §2.3). 

The resulting H-band light curve is included in 
Fig. 1, and a comparison of all three IR band light 
curves (J, H, K/K') is shown in Fig. 5. Generally, 
we find moderate variability of A(J, H, K/K 1 ) < 
0.2 within ~ a few days. The various infrared 
bands trace each other very well, and there does 
not appear to be any significant time lags be- 
tween the different bands. It is obvious that there 
are issues with the cross calibration between the 
three observatories. In particular, the Mt. Abu 
observations (using the J, H, and K' filters) in 
mid-December of 2003 indicate a significant drop 
of the average level of near-infrared flux. While 
there appears to be a similar flux drop around the 
same time in the radio light curves, the much more 
densely sampled optical light curves do not indi- 



cate a similar feature. However, it is obvious that 
the near-IR data around this time exhibit an un- 
usually large scatter (by ^0.5 mag on sub-hour 
time scales), they may be affected by calibration 
uncertainties. We note that rebinning of the IR 
data in time bins of up to 30 minutes did not im- 
prove the quality of the data because of the ex- 
treme nature of the apparent flux fluctuations. We 
also point out that the use of the K' filter at Mt. 
Abu vs. the K filter at the other two IR observa- 
tories, could only explain offsets of less than the 
individual measurement errors, so this effect was 
neglected in our analysis. 

2.3. Optical observations 

In the optical component of the extensive 
WEBT campaign, 24 observatories in 15 coun- 
tries contributed 7611 individual photometric data 
points. The observing strategy and data analysis 
followed to a large extent the standard procedure 
described previously for a similar campaign on 
BL Lacertae in 2000 (Villata et al. 2002). For 
more information about WEBT campaigns see 
also: Villata et al. (2000, 2002, 2004a,b); Raiteri 
ct al. (2001, 2005) 

It had been suggested that, optimally ob- 
servers perform photometric observations alter- 
nately in the B and R bands, and include com- 
plete (U)BVRI sequences at the beginning and the 
end of each observing run. Exposure times should 
be chosen to obtain an optimal compromise be- 
tween high precision (instrumental errors less than 
~ 0.03 mag for small telescopes and ~ 0.01 mag 
for larger ones) and high time resolution. If this 
precision requirement leads to gaps of 15 - 20 min- 
utes in each light curve, we suggested to carry out 
observations in the R band only. Observers were 
asked to perform bias (dark) corrections as well 
as flat-fielding on their frames, and obtain instru- 
mental magnitudes, applying either aperture pho- 
tometry (using IRAF or CCDPHOT) or Gaussian 
fitting for the source 3C 66A and the compari- 
son stars no. 13, 14, 21, and 23 in the tables of 
Gonzalez-Perez et al. (2001), where high-precision 
standard magnitudes for these stars have been 
published. This calibration has then been used 
to convert instrumental to standard photometric 
magnitudes for each data set. In the next step, un- 
reliable data points (with large error bars at times 
when higher-quality data points were available) 
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were discarded. Our data did not provide evi- 
dence for significant variability on sub-hour time 
scales. Consequently, error bars on individual 
data sets could be further reduced by re-binning 
on time scales of typically 15 - 20 min. 

Finally, there may be systematic offsets be- 
tween different instruments and telescopes. Wher- 
ever our data sets contained sufficient indepen- 
dent measurements to clearly identify such off- 
sets, individual data sets were corrected by ap- 
plying appropriate correction factors. In the case 
of BL Lacertae and similar sources (e.g., Villata 
et al. 2002), such corrections need to be done on 
a night-by-night basis since changes in the seeing 
conditions affect the point spread function and 
thus the precise amount of the host galaxy con- 
tamination. However, in the absence of a signif- 
icant host galaxy contribution, these systematic 
effects should be purely instrumental in nature 
and should thus not depend on daily seeing con- 
ditions. Thus, we opted to introduce only global 
correction factors for entire single-instrument data 
sets on 3C 66A. This provided satisfactory results 
without obvious residual inter-instrumental incon- 
sistencies. The resulting offsets are listed in Tab. 
4. 

In order to provide information on the intrin- 
sic broadband spectral shape (and, in particular, 
a reliable extraction of B - R color indices), the 
data were de-reddened using the Galactic Extinc- 
tion coefficients of Schlegel et al. (1998), based on 
A B = 0.363 mag and E(B - V) = 0.084 mag 4 . 
As mentioned in the introduction, the R magni- 
tude of the host galaxy of 3C 66A is ~ 19 mag, so 
its contribution is negligible, and no host-galaxy 
correction was applied. 

2.3.1. Light curves 

As a consequence of the observing strategy de- 
scribed above, the R- and B-band light curves are 
the most densely sampled ones, resulting in sev- 
eral well-sampled light curve segments over ~ 5 - 
10 days each, with no major interruptions, except 
for a gap of a few hours due to the lack of cov- 
erage when 3C 66A would have been optimally 
observable from locations in the Pacific Ocean. 
The R-band light curve over the entire duration 
of the campaign is included in Fig. 1 and com- 
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Fig. 5. — Light curves in the near-infrared J, H, 
and K bands over the entire duration of the cam- 
paign. 
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Fig. 6. — Light curves in the optical U, B, V, R, 
and I bands over the entire duration of the cam- 
paign. 
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Table 4 

Inter-instrumental correction offsets for optical data 



Observatory 


B 


V 


R 


1 


Armenzano 


-\-.02 




+ .02 




Bell 










Boltwood 


no 








Catania 










Crimea 


_ 


+ .03 


_ 


_ 


Heidelberg 










MDM 










Michael Adrian 










Mt. Lemmon 










Mt. Maidanak 


+ .01 








Nyrola 






-.03 




Osaka 










Perugia 




+ .08 


+ .05 


+ .04 


Roque(NOT) 










Sabadell 










San Pedro Martir 


+ .08 




+ .07 




Shanghai 






-.03 




Skinakas 


-.01 


-.04 


-.04 


-.05 


Sobacksan 










St. Louis 


-.06 








Torino 


-.09 






-.03 


Tuorla 




+ .02 


-.02 





pared to the light curves at all other optical bands 
in Fig. 6. These figures illustrate that the object 
underwent a gradual brightening throughout the 
period July 2003 - February 2004, reaching a max- 
imum at R w 13.4 on Feb. 18, 2004, followed by a 
sharp decline by AR ~ 0.4 mag within ~ 15 days. 
On top of this overall brightening trend, several 
major outbursts by AR <~ 0.3 - 0.5 mag on time 
scales of <~ 10 days occurred. The two most dra- 
matic ones of these outbursts peaked on July 18, 
2003, and Nov. 1, 2003. Details of the Nov. 1 
outburst are displayed in Figs. 7-9, respec- 
tively. We find evidence for intraday microvari- 
ability of AR ~ 0.05 mag on time scales down to 
~ 2 hr. One example for such evidence is illus- 
trated in Fig. 10 which shows the intranight data 
from the Torino Observatory for JD 2,452,955 (= 
Nov. 11, 2003). The figure presents the instru- 
mental magnitudes of 3C 66A, comparison star 
A, and the difference of both. While the differ- 
ence seems to be affected by seeing effects at the 
beginning and end of the observation (rising and 
falling in tandem with the instrumental magni- 
tudes), there is evidence for intraday variability 
around ~ JD 2,452,955.5, where there were only 
very minor changes in the atmospheric opacity. 

Visual inspection of individual major outbursts 
suggests periods of more rapid decline than rise. 
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Fig. 7. — Details of the R-band light curve during 
the major outburst around Nov. 1, 2003. The 
boxes indicate the light curve segments shown in 
more detail in Figs. 8-9. 
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Fig. 8. — Details of the R-band light curve during 
the rising phase of the major outburst around Nov. 
1, 2003. 



2.3.2. Periodicity analysis 

For several blazars, including 3C 66A, period- 
icities on various time scales, from several tens 
of days, to several years, have been claimed (see, 
e.g. Rieger 2004, and references therein for a more 
complete discussion) . In order to search for a pos- 
sible periodicity in the optical R-band light curve 
from our WEBT campaign, we have performed a 
Fourier analysis of the R-band light curve after 
interpolation of the light curve over the (short) 
gaps in the available data. Such an analysis be- 
comes unreliable for periods of more than ~ 1/10 
of the length of the data segment. Consequently, 
our data can be used for a meaningful analysis on 
time scales of r < 30 - 50 d. No evidence for 
a periodicity in this range has been found. How- 
ever, we note a sequence of large outbursts on July 
18 (MJD 52838), Sept. 5 (MJD 52887), Nov. 1 
(MJD 52944), Dec. 28 (MJD 53001), 2003, and 
Feb. 18, 2004 (MJD 53053), which are separated 
by intervals of 49, 57, 57, and 52 days, respec- 
tively. These intervals appear remarkably regular, 
and only slightly shorter than the 65-day period 
claimed by Lainela et al. (1999). The duration 
of our monitoring campaign is too short to do a 
meaningful analysis of the statistical significance 
of this quasi-periodicity. This question will be re- 
visited in future work, on the basis of a larger data 
sample, including archival optical data. 
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Fig. 9. — Details of the R-band light curve during 
the decaying phase of the major outburst around 
Nov. 1, 2003. 



2.4. X-ray observations 

3C 66A was observed by the Rossi X-Ray Tim- 
ing Explorer (RXTE) 26 times between Septem- 
ber 19 and December 27, 2003, for a total ob- 
servation time of approximately 200 ks. Analy- 
sis of the RXTE PC A data was carried out using 
faint-source procedures as described in the RXTE 
Cookbook. Standard selection criteria were used 
to remove disruptions from the South Atlantic 
Anomaly, bright Earth, instrument mispointing, 
and electron contamination. For the entire data 
set, only PCUs and 2 were activated, but data 
from PCU was discarded due to its missing 
propane veto. There is also a persistent problem in 
the background model of the PCA around the 4.78 
keV Xenon L-cdgc, which causes a small anomaly 
in count rates near this energy. For faint sources 
like 3C 66 A, this anomaly can be significant, so 
photons in the three energy bins surrounding this 
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energy were discarded. 

Photons from 3-10 keV were extracted from 
the data using only the top layer of PCU2. Ex- 
tracting higher-energy photons proved unreliable, 
as the source flux was close to the PC A source- 
confusion limit. A spectrum was compiled using 
data from all 26 source observations. The spec- 
trum was fit to a simple power-law, which resulted 
in an energy spectral index of a = 1.47 ± 0.56 
and a flux $i koV = (4.06 ± 0.49) x 10~ 3 pho- 
tons cm~ 2 s _1 keV -1 , with x 2 / n = 13.6/14. 
Galactic absorption was not included in the fit 
model, as the lack of data below 3 keV precludes 
a useful constraint on N# . A plot of the spectrum 
and fit is shown in Figure 11. 

The PCA is a non-imaging instrument with a 
FWHM response of 1°, so there is the potential 
that the signal found here is a sum over multiple 
sources. Most notably, the FR-1 galaxy 3C 66B 
lies only 6' from 3C 66A and is a known X-ray 
emitter. Though indistinguishable from 3C 66A 
with the PCA, its typical flux is at a level of ~ 2 % 
of the measured PCA flux of 3C 66A (Croston et 
al. 2003), and is thus not a significant contribu- 
tor. In the 3EG catalog, there is also confusion of 
3C 66A with a nearby pulsar, PSR J0218+4232. 
At a source separation of 58', the PCA response 
for the pulsar is low for 3C 66A pointed obser- 
vations, and thus it does not make a significant 
contribution. 

The 3-10 keV photons were binned into 24- 
hour periods to render a light curve, which is in- 
cluded in Fig. 1. No significant flux variations 
or transient events were seen. To check for spec- 
tral variations, the data set was divided into four 
roughly equal sub-periods, and power-law spectra 
fitted to each sub-period. Fits to all four sub- 
periods were consistent with the spectral parame- 
ters obtained from the entire data set, and all sub- 
periods were consistent with each other. From this 
we conclude that there were no variations of cither 
flux or spectral shape in the X-ray band during the 
observing period within the detection capabilities 
of RXTE. 

2.5. Gamma-ray observations 

At VHE 7-rays, 3C 66A was observed con- 
temporaneously with our broadband campaign by 
the Solar Tower Air Cerenkov Effect Experiment 
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shows clear evidence for intranight variability 
during the central portion of this observation 
(~ JD 2,452,955.4 - - JD 2,452,955.5). 
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Fig. 11 — Average RXTE spectrum (3 - 10 keV) 
and power-law fit for all observations of 3C 66A 
taken as part of the multiwavelength campaign. 
See text for best-fit parameters. 
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(STACEE) and by the 10-m Whipple Telescope of 
the Very Energetic Radiation Imaging Telescope 
Array System (VERITAS) collaboration. 

STACEE took a total of 85 28-minutc on-off 
pairs of data, totalling 33.7 hours of live time on- 
source. After data-quality cuts and padding, a 
total of 16.3 hours of on-source live time remained. 
A net on-source excess of 1134 events attributed to 
photons of energy E > 100 GeV was seen against 
a background of 231,742 events. The details of the 
STACEE data analysis are published in a separate 
paper (Bramel et al. 2005). 

The source excess quoted above corresponds to 
a significance of 2.2 a, which is insufficient to claim 
a detection, but can be used to establish firm up- 
per flux limits. For various assumptions of an 
underlying intrinsic source power-law spectrum, 
those upper limits are listed in Table 5. The data 
were also binned in 24-hr segments. This revealed 
no significant excess on any given day, and thus 
no evidence for a statistically significant transient 
event during the campaign period. 

The VERITAS collaboration observed 3C 66A 
during the period Sept. 27, 2003 - Jan. 13, 2004, 
obtaining a total of 31 on-off pairs with typically 
~ 27.6 min. of on-source live time per pair. The 
data were analyzed following the standard Whip- 
ple data analysis procedure (see, e.g. Falcone et 
al. 2004). The log of the individual observations 
and 99.9 % confidence upper limits, assuming an 
underlying Crab-like source spectrum with an un- 
derlying energy index of 2.49 (Hillas et al. 1998), 
are listed in Table 6. Combining all measure- 
ments results in a 99.9 % confidence upper limit 
of 0.91 x 10~ n ergs cm~ 2 s -1 . 

3. Spectral variability 

In this section, we will describe spectral vari- 
ability phenomena, i.e. the variability of spec- 
tral (and color) indices and their correlations with 
monochromatic source fluxes. As already men- 
tioned in the previous section, no evidence for 
spectral variability in the X-ray regime was found. 
Consequently, we will concentrate here on the op- 
tical spectral variability as indicated by a change 
of the optical color. In particular, our observing 
strategy was optimized to obtain a good sampling 
of the B - R color index as a function of time. 

Based on our observation (see §2.3) of no sig- 
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Fig. 12. — Light curve of the B and R magni- 
tudes and the B - R color index of 3C 66A over 
the duration of the entire campaign. While the B 
and R light curves are very well correlated, with 
no significant detection of a time delay, maxima 
of the spectral hardness (minima of B - R) sys- 
tematically precede B and R outbursts by a few 
days. 



0.6 



0.7 



; o.s 



0.9 




+ 



14.5 



_l I l_ 



14 

R magnitude 



13.5 



Fig. 13. — Optical hardness - intensity diagram 
for the complete data set over the entire dura- 
tion of the campaign. While a positive hardness- 
intensity correlation seems to be present at low 
brightness, R > 14.0, no clear correlation can be 
found in brighter states. 
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Table 5 

99 % confidence upper limits on the vhe photon flux from the stacee observations for 
various assumptions of an underlying intrinsic photon spectral index t = a + 1 



r 


Sthr [GoV] 


dN/dE at B thr [photons m~ 2 s -1 GoV -1 ] 


2.0 


200 


5.23 x 10~ 9 


2.5 


184 


9.39 x 10~ 9 


3.0 


150 


2.26 x 10~ 8 


3.5 


147 


3.10 x 10~ 8 



Table 6 

Observation log and 99 % confidence upper limits on the integrated > 390 GeV photon 
flux from the whipple observations of 3c 66a, assuming a crab-like spectrum. 



Date 


Start time [UTC] 


On-Source Live Time [min.] 


99.9 % UL [lO" 11 ergs cm- 2 s" 1 ] 


030927 


10:57 


27.73 


2.05 


030928 


09:04 


27.70 


3.16 


030929 


08:16 


27.75 


2.73 


031002 


08:58 


27.66 


1.77 


031002 


09:27 


27.73 


2.35 


031002 


10:25 


27.69 


1.82 


031023 


06:58 


27.55 


6.19 


031024 


07:07 


27.57 


5.65 


031028 


06:40 


27.55 


6.35 


031028 


07:08 


27.57 


7.00 


031029 


06:41 


27.57 


6.62 


031029 


08:15 


19.99 


5.29 


031030 


06:27 


27.61 


4.39 


031101 


07:17 


27.58 


4.51 


031125 


05:49 


27.55 


6.32 


031128 


04:42 


27.57 


2.81 


031214 


03:06 


27.62 


4.87 


031217 


03:00 


27.63 


5.95 


031217 


04:08 


27.54 


5.23 


031217 


04:37 


27.58 


5.23 


031217 


05:05 


27.55 


6.71 


031217 


05:33 


27.58 


3.67 


031220 


03:38 


27.63 


2.60 


031220 


05:09 


27.60 


5.95 


031220 


06:07 


27.66 


4.45 


031223 


02:48 


27.70 


4.63 


031223 


04:01 


27.71 


4.03 


031223 


04:31 


27.68 


3.88 


031223 


04:59 


27.70 


2.94 


031225 


04:35 


27.64 


5.44 


040113 


02:17 


27.58 


5.23 
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nificant variability on time scale less than ~ 1 hr, 
we extracted B - R color indices on measurements 
of B and R magnitudes separated by less than 20 
min. of each other. Fig. 12 shows a comparison 
of the light curves of the B and R magnitudes and 
the B - R color index. The figure reveals an inter- 
esting new result: While the B and R light curves 
are well correlated, with no significant detection 
of a time lag from a discrete correlation function 
analysis (see §4) , maxima of the spectral hardness 
(i.e., minima of the value of B - R) seem to pre- 
cede B and R band flux outbursts systematically 
by a few days. In other words, the optical spec- 
tra harden at the onset of a major outburst, and 
continue to soften as the flare progresses through 
both the rising and the decaying phase. However, 
Fig. 12 also indicates that the apparent lead time 
of the B - R hardness maximum varies within a 
range of a few days. This might be the reason that 
a cross-correlation analysis between the B - R in- 
dex and the B or R light curves does not reveal a 
strong peak. In fact, the DCF between those time 
series does never exceed values of ~ 0.3 at all. 

Very interesting and intriguing is also the B - 
R behaviour of the source around JD 2452955 - 
JD 2452995, where no major flux outbursts are 
observed, but several episodes of significant spec- 
tral hardening by A(B — R) > 0.1 (corresponding 
to spectral-index variations of Aa op t > 0.2) are 
observed. 

In Fig. 13 we have plotted the B - R index vs. 
R-band magnitude for the entire data set. The 
plot shows that there is a weak indication of a 
positive hardness-intensity correlation at low flux 
states with R > 14.0. At higher flux levels, no 
correlation is apparent, which might be a conse- 
quence of the result found above, that the B - R 
hardness actually peaks during the rising phase of 
individual outbursts. This is in contrast with re- 
cent results of Vagnetti et al. (2003) who, on the 
basis of a smaller data set, have found a consistent 
trend of B - R hardening with increasing B-band 
flux, independent of the actual flux value. 

The optical spectral hysteresis visible in our 
data is further illustrated in Fig. 14, where we 
have plotted the B - R vs. R hardness- intensity di- 
agram for the example of one individual outburst 
on Nov. 1, 2003 (MJD 52944). The top panel 
shows the hardness-intensity diagram for the en- 
tire flare, which is then split up into the rising and 



the decaying phase of the flare in the middle and 
lower panels, respectively. The figure illustrates 
that as the flare rises up to the peak R bright- 
ness, the spectrum already begins to soften sig- 
nificantly as R < 13.6. Similar trends are found 
for other optical flares as well. Possible physical 
interpretations of this trend, along with detailed 
modelling of the SED and spectral variability will 
be presented in a separate paper ( Joshi & Bottchcr 
2005, in preparation). 

4. Inter-band cross-correlations and time 
lags 

In this section, we investigate cross-correlations 
between the measured light curves at different 
frequencies, within individual frequency bands as 
well as broadband correlations between different 
frequency bands. Because of relatively poor sam- 
pling of the radio and IR light curves, our searches 
for correlations between different radio and IR 
bands and between those wavelength bands with 
the optical ones on short time scales compared 
to the duration of our campaign did not return 
conclusive results. Also, in agreement with our 
non-detection of a periodicity on time scales of 
< 50 days, a discrete autocorrelation function 
(DACF Edelson & Krolik 1988) analysis of indi- 
vidual optical light curves did not return results 
beyond artificial "periodicities". Such artifacts 
can be attributed to the quasi-periodic, uneven 
time coverage which results in spurious detections 
of periodicities at multiples of days and multiples 
of ~ 0.5 d (the time delay between the peak cover- 
age from the heavily contributing observatories in 
Europe and in East Asia). A discrete autocorre- 
lation function analysis of the R-band light curve 
also revealed secondary peaks at ±4 days with a 
correlation coefficient close to 1. 

The cleanest results of a DCF analysis are obvi- 
ously expected for the most densely sampled light 
curves, which we obtained in the B and R bands. 
We performed DCFs between the R and B band 
light curves on a variety of time scales with a va- 
riety of binning intervals, ranging from 15 min. to 
10 days. On intraday time scales, we find consis- 
tently a sharp peak at delay with all bin sizes 
we used, which indicates no evidence for a time 
delay between variability patterns in the R and 
B band reaching their peak fluxes in each band. 
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The DCFs on short time scales are dominated by 
the artificial 0.5 d and 1 d periodicities mentioned 
above. The discrete correlation function between 
R and B band magnitudes on time scales of sev- 
eral days seems to indicate a strong correlation at 
t ~ — 4 days, corresponding to a lead of the B vs. 
the R band. However, this might be a consequence 
of the probably artificial 4-day "periodicity" of the 
R-band light curve discussed in the previous para- 
graph, which prevents us from making any claim 
about the detection of a 4-day delay between the 
R and B band light curves on the basis of our DCF 
analysis. 

5. Broad-band spectral energy distribu- 
tions 

From the various flux measurements described 
in detail in the previous sections, we can now com- 
pose contemporaneous spectral energy distribu- 
tions (SEDs) of 3C 66A at various times during 
our campaign. Significant variability was only de- 
tected in the radio, IR, and optical bands. In those 
bands, we extracted SEDs for four epochs: Dur- 
ing two major outbursts, around Nov. 1, 2003, 
and Dec. 28, 2003, a minor outburst around Oct. 
1, 2003, as well rather quiescent state around 
Nov. 11, 2003. De-reddened optical and near- 
IR magnitudes were converted to fluxes using the 
zero-point normalizations of Bessel, Castelli, & 
Plez (1998). The resulting SEDs are plotted in 
Fig. 15. Due to the relatively poor sampling, ra- 
dio and IR data were often not quite simultane- 
ous with the optical spectra, which were extracted 
near the peaks of the invididual outbursts men- 
tioned above. In this case, the closest near-IR and 
radio data points were chosen. This led to time off- 
sets between the optical and the radio data points 
of up to ~ 5 d. Given the relatively long time 
scale and moderate amplitude of variability at ra- 
dio wavelengths, we are confident that this did not 
introduce serious distortions of the low-frequency 
SED. 

In addition to the data taken during our 
2003/2004 campaign, we have included histori- 
cal X-ray measurements, to indicate the degree 
of X-ray variability observed in this source and 
the historical average GeV 7-ray flux measured 
by the EGRET instrument on-board the Comp- 
ton Gamma-Ray Observatory from 5 observations 



between Nov. 1991 and Sept. 1995 (Hartman et 
al. 1999). 

The shape of the time-averaged optical (U)BVRI 
spectra, together with the very steep X-ray spec- 
tral index, indicates that the vF v peak of the 
synchrotron component of 3C 66A is typically 
located in the optical range. The shape of our 
best-fit RXTE spectrum provides strong evidence 
that the synchrotron component extends far into 
the X-ray regime and intersects the high-energy 
component in the time-averaged SED of 3C 66A 
at > 10 keV during our core campaign period. 
Due to the lack of simultaneous GeV 7-ray cov- 
erage and of a firm detection at > 100 GeV, we 
can not make a precise statement concerning the 
level of 7-ray emission during our campaign. How- 
ever, if the historical EGRET flux is representative 
also for the time of our campaign, then the total 
energy output in the low-frequency and the high- 
frequency components of the SED of 3C 66A are 
comparable, as is typical for the class of interme- 
diate and low-frequency peaked BL Lac objects. 

6. Generic parameter estimates 

In this section, we discuss some general con- 
straints on source parameters. We will first (§6.1) 
focus on parameter estimates using the spectral 
energy distribution and optical intraday variabil- 
ity measurements, most relevant to the innermost 
portion of the jet outflow, closest to the central en- 
gine. In §6.2, we will use the results of our VLBA 
observations to estimate parameters of the rela- 
tivistic flow on parsec scales. 

6.1. Parameters of the inner jet 

From the minimum variability time scale of 
Ai m i n ~ 2 hr, we can estimate the size of the 
emitting region as R < cDAt m i n , where D = 
(r [1 — (3r cos 6* bs]) _1 , where T is the bulk Lorentz 
factor of the emitting region, /3r c is the cor- 
responding speed, and o bs is the observing an- 
gle. This yields R < 2.2 x 10 15 Si cm, where 
S = 10 (5i ~ 15 as an estimate from the limits on 
superluminal motion and from previous modeling 
efforts, as mentioned in the introduction. 

An estimate of the co-moving magnetic field 
can be found by assuming that the dominant por- 
tion of the time-averaged synchrotron spectrum 
is emitted by a quasi-equilibrium power-law spec- 
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Fig. 14. — Optical hardness - intensity diagram 
for the flare around Nov. 1, 2003. The opti- 
cal spectrum already begins to soften significantly 
while the R-band flux is still rising at R < 13.6. 
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trum of electrons with N e (j) — uq Vb 7~ p for 
7i < 7 < 72! here, Vb is the co- moving blob 
volume. Based on the X-ray spectral index of 
a « 1.5, we find a particle spectral index of p w 4. 
Since the synchrotron cooling time scale for X- 
ray emitting electrons might be much shorter than 
the dynamical time scale (see Eq. 3), the X-rays 
are likely to be produced by a cooled electron dis- 
tribution. In this case, the index p m 4 corre- 
sponds to a distribution of electrons injected into 
the emission region with an original index q « 3. 
The normalization constant no is related to the 
magnetic field through an equipartition parame- 
ter es = us/u e (in the co- moving frame). Note 
that this equipartition parameter only refers to 
the energy density of the electrons, not account- 
ing for a (possibly dominant) energy content of 
a hadronic matter component in the jet. Un- 
der these assumptions, the magnetic field can be 
estimated as described, e.g., in Bottcher et al. 
(2003). Taking the vF v peak synchrotron flux /| y 
at the dimensionless synchrotron peak energy e sy 
as ~ 5 x 10~ n ergs cm" 2 s" 1 at e sy « 5 x 10~ 6 
and R w 3.3 x 10 15 cm, we find 



B, 



4.4 5- 



-12/7 
e 73 



G. 



(1) 



2.9e 2 R /7 G for 5 = 15. 



which yields B e 

We can further use this magnetic-field value to 
estimate the end points of the electron spectrum 
since the low-energy end, 71 might correspond to 
the vF v peak of the synchrotron spectrum, and 
the synchrotron high-energy cutoff corresponds to 
72. Generally, we find 



Fig. 15. — Time-averaged spectral energy distri- 
butions of 3C 66A at various epochs during the 
core campaign period in 2003. The turquoise di- 
amonds indicate the average of the RATAN-600 
measurements from Oct. 11 and 14, 2003. His- 
torical soft X-ray spectra are plotted for compar- 
ison; the hard X-ray spectrum labelled "RXTE 
2003" indicates the result of a power-law fit to the 
RXTE PCA observations of this campaign. The 
historical (non-contemporaneous) average of the 5 
EGRET pointings on the source is also included 
for reference. All data that are not contempora- 
neous with our campaign are indicated by dotted 
curves. 



where ^15 is the characteristic synchrotron fre- 
quency in units of 10 15 Hz. With our standard 
parameters, this yields 71 rj 3.1 x 10 3 , and 72 = 
1.5 x 10 5 , if the synchrotron cutoff occurs around 
10 keV. We can also use this to estimate the syn- 
chrotron cooling time scale of electrons in the ob- 
server's frame: 

(3) 
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For optical frequencies, this yields observed cool- 
ing time scales of the order of ~ 2 hr, in agreement 
with the observed minimum variability time scale. 
This, however, raises an important caveat: The 
observed minimum variability time scale may, in 
fact, be a reflection of the electron cooling time 
scale rather than the dynamical time scale, as we 
had assumed when choosing our estimate for the 
source size R. For this reason, a more detailed fu- 
ture investigation of possible short-term variabil- 
ity at X-ray frequencies will be extremely impor- 
tant to resolve this issue. If X-ray variability on 
shorter time scales than ~ 2 hr is found, the emis- 
sion region would then have to be more compact 
than the R <~ 3 x 10 15 cm which we had estimated 
here, and X-ray spectral hysteresis patterns could 
arise. In contrast, if the X-ray variability time 
scale is found to be consistent with the optical 
one, then it would have to be dominated by the 
dynamical time scale and thus be largely achro- 
matic. Thus, in that case no significant X-ray 
spectral hysteresis would be expected. 

6.2. Parameters of the pc-scale outflow 

In our analysis of the 22 GHz and 43 GHz 
VLBA maps of 3C 66A, we had found a rather 
smooth jet, with only one of the 6 Gaussian com- 
ponents (CI) showing evidence for superluminal 
motion of /3 app = (8.5 ± 5.6) hr 1 . If we consider 
the measured speed of the component CI to be 
close to the maximum as observed under the su- 
perluminal angle given by cos(6^ s ) = 0r, we can 
estimate the minimum Lorentz factor and maxi- 
mum angle between the jet and our line of sight. 
With h = 0.7 and by taking the lower limit of the 
/?app > (8-5 — 5.6) h^ 1 = 2.9 h^ 1 , we have a lower 
limit to the bulk Lorentz factor of T > 4.3 and an 
upper limit to the viewing angle of ^obs < 27.2°. 
The large error in /3 app results in very weak con- 
straints on T and 6> bs, and they should be consid- 
ered as very conservative estimates. Much more 
accurate values should be available after the anal- 
ysis of all 9 epochs of the VLBA monitoring pro- 
gram (Savolainen et al. 2005, in preparation). 

From the measured fluxes at 22 and 43 GHz, 
we have calculated the brightness temperatures 
along the jet. The largest brightness tempera- 
tures are found for the core, yielding lower limits 
of T b COTC ~ 10 11 K (see Tab. 2). Here, we have 
to mention the caveat that the errors for the size 



10" 



10' 





\^ s 22 - -2.2HJ.2 


















ft 










■ 



























10"' 10° 10' 

Distance from the core [mas] 

Fig. 16. — VLBA brightness temperatures as a 
function of distance from the core, at 22 GHz 
(filled symbols) and 43 GHz (open symbols), 
for the three epochs: 2003.78 (circles), 2003.83 
(squares), 2004.08 (stars). The solid lines are our 
best-fit power-laws to the brightness-temperature 
profiles; the vertical dashed lines indicate the jet 
regions B2 and B3, where a significant deviation 
from the power-law profile is found, coincident 
with an apparent kink in the jet flow direction. 
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of the core are very difficult to estimate, since our 
model-fitting yields very asymmetric probability 
density distributions for the core size. Essentially, 
the data is also consistent with a point source core, 
so we can only give an upper limit to the core 
size. Thus, no upper limit to the core brightness 
temperature is available. Our calculated best-fit 
core brightness temperatures and corresponding- 
lower limits are listed in Tab. 2 and can be com- 
pared with an equipartion brightness temperature 
(Readhead 1994) for 3C 66A of T b , cp = 6 * 10 10 K. 
If the equipartition assumption holds, the mini- 
mum Dopplcr factor is around 2. The same result 
was obtained by Lahteenmaki & Valtaoja (1999) 
who used variability arguments to calculate the 
Doppler factor. Because we do not have upper 
limits for the brightness temperature, we cannot 
calculate a more precise value for the Doppler fac- 
tor. However, if the core brightness temperatures 
calculated from the best model-fit values arc not 
too far from the truth, the core Doppler factor 
probably lies somewhere between 2 and 15 with 
the average best-fit value of 6. 

We have also calculated the brightness temper- 
atures for components other than the core at both 
22 and 43 GHz. In Fig. 16, the component bright- 
ness temperatures are plotted against the compo- 
nent distance from the core. The figure illustrates 
that the points form two well-defined power laws 
with a different scaling at different frequencies. All 
three epochs are included in the plot (and coded 
by colors). The brightness temperature gradient 
along the jet is well described by a power law 
Tb oc r s with s = — 2. There is an exception to 
the power law between 1.8 and 2.3 mas from the 
core, where the brightness temperature increases 
above the power law fit by a factor of 3 - 5. This 
takes place in the jet region B, where the jet is 
clearly bending. This brightening could be a re- 
sult of a temporary increase of the Doppler factor 
caused by the bending of the jet towards our line 
of sight. An alternative interpretation could be 
given by impulsive particle acceleration in the re- 
gion B by a standing shock wave caused by the 
interaction of the continuous, relativistic outflow 
with the observed kink. 

The power law gradient of the brightness tem- 
perature along the jet suggests that 3C 66A ex- 
hibits a smooth pc-scale jet without any promi- 
nent particle acceleration sites (shocks) other than 



the core and possibly the bend in region B. If 
we assume that the magnetic field B, the elec- 
tron density N and the cross-sectional diameter 
D of the jet can also be described by power laws: 
N (x r n , B oc r b and D oc r d , the brightness tem- 
perature is expected to fall with r as Jf, oc r s (with 
s < 0). If optically thin synchrotron emission and 
a constant Lorentz factor of the emitting electrons 
are assumed, s = d + n + 6(1 + a), where a is the 
spectral index (see, e.g. Kadler et al. 2004). If we 
know the values of s, d and a, we can calculate 
the relationship between b and n, or, by assuming 
equipartition, the actual values of b and n. We 
know now that s — — 2, and the average spectral 
index is (022-43) = 0.15. 

We still need to find out the value of d. d = 1 
would correspond to a constant opening angle 
(conical jet), but our imaging results allow us to 
estimate a more precise value from the data by 
plotting the component size versus distance from 
the core (see Fig. 17). Surprisingly, up to 2.5 mas 
from the core, the jet is very tightly collimated. 
A power-law fit to the observed correlation be- 
tween the lateral jet size and distance from the 
core yields d — 0.6. After 2.5 mas, the jet opens 
up (or disrupts) and d becomes larger than 1. 
This might be another indication that the angle 
between the jet and our line of sight is changing 
along the jet. The observed opening angle of the 
jet, V'obs, is related to the intrinsic opening angle 
V>int, through i/'obs ~ V'int/sin^obs) for small val- 
ues of ?/>i n t and 6> bs- Thus, the observed opening 
angle increases when the jet bends towards our line 
of sight, which is in accordance with the observed 
brightening of the jet in the region B. 

If equipartition holds, we may assume n = 2b. 
Then, taking d = 0.6 yields b = —0.8 and n = 
— 1.7. For the case of a conical jet (d — 1), we find 
n = —1.9 and b = —1.0. If the equipartition as- 
sumption holds, these results mean that in a con- 
ical jet, the magnetic field decays as B oc D -1 , 
and in a collimated jet [d — 0.6), it decays as 
B oc D~ 13 , implying a magnetic field that is pre- 
dominantly transverse to the jet axis. This is in 
accordance with other observational results on the 
magnetic-field orientation in BL Lac objects in 
general (see, e.g. Gabuzda et al. 2004). A mag- 
netic field being predominantly transverse to the 
jet agrees well with overall toroidal field configura- 
tion which is expected in Poynting flux dominated 
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(PFD) jet models. Another possibility for creat- 
ing predominantly transverse magnetic field is the 
enhancement of the transverse magnetic field com- 
ponent by a series of shocks, but this seems to be 
less likely in this case, since we do not observe 
bright knots in the jet, but a rather smooth flow. 

7. Summary 

We have observed 3C 66A in a massive multi- 
wavelength monitoring campaign from July 2003 
to April 2004. Monitoring observations were car- 
ried out at radio, infrared, optical, X-ray, and 
VHE 7-ray observations. The main observational 
results of our campaign are: 

• In the optical, several outbursts by Am ~ 
0.3-0.5 over time scales of several days were 
observed. 

• Optical intraday microvariability (Am < 
0.05) on time scales of Ai var ~ 2 hr was 
detected. 

• No clear evidence for periodicity was found, 
but a quasi-regular sequence of several ma- 
jor outbursts separated by ~ 50 - 57 d was 
identified. 

• Large optical flares (on time scales of sev- 
eral days) seem to exhibit optical spectral 
hysteresis, with the B - R hardness peaking 
several days prior to the R- and B-band flux 
peaks. 

• The 3-10 keV X-ray spectrum is best fitted 
with a single power-law with energy index 
a = 1.47 ± 0.56, indicating that the transi- 
tion between the synchrotron and the high- 
energy components occurs at photon ener- 
gies of > 10 keV. 

• Radio VLBA monitoring observations reveal 
a rather smooth jet with no clearly dis- 
cernible knots or hot spots except for mod- 
erate brightening in the region B where the 
jet is bending. 

• Decomposition of the VLBA radio structure 
into Gaussian components reveals superlu- 
minal motion only for one out of 6 compo- 
nents. Its apparent speed is /3 app = (8.5 ± 
5.6) h- 1 = (12.1 ± 8.0) for h = 0.7. 



• The radio brightness temperature profile 
along the jet, along with its observed geome- 
try, D oc r 0,6 , suggest a magnetic-field decay 
B oc Z?~ 13 , indicating a predominantly per- 
pendicular magnetic-field orientation. 

• STACEE observations revealed a 2.2 a ex- 
cess, which provided strict upper limits at 
E-y > 150 GcV. Additional VHE 7-ray limits 
at E 1 > 390 GeV resulted from simultane- 
ous Whipple observations. 

• The broadband SED of 3C 66A during our 
campaign suggests that the synchrotron 
component peaks in the optical and extends 
far into the X-ray regime, out to at least 
- 10 keV. 

• The following parameters of the synchrotron 
emission region near the core can be esti- 
mated: 

5 « 15 

R w 3.3 x 10 15 cm 
B « 2.9e 2 J 7 G 

71 w 3.1 x 10 3 

72 w 1.5 x 10 5 

p « 4 (4) 

• No X-ray variability was detectable by 
RXTE on a ~ 1 day time scale. Future 
X-ray observations with more sensitive X- 
ray detectors will be important to probe for 
rapid X-ray variability and X-ray spectral 
hysteresis in order to put more stringent 
constraints on the source size, the nature of 
the variability mechanism, and the compo- 
sition and energetics of the emitting plasma 
in the jet. 
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Fig. 17. — Transverse diameter D (FWHM of 
Gaussian component fits) of the jet as a function 
of distance from the core for the inner 3 mas. The 
symbols are the same as in Fig. 16. The solid 
curve illustrates the best-fit power law D oc r d , 
which yields d — 0.6. 
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